Electron-phonon interaction in cuprate-oxide superconductors 
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We propose a novel electron-phonon interaction arising from the modulation of the superexchangc 
interaction by phonons. It is enhanced by spin and superconducting fluctuations, which are devel- 
oped mainly because of the superexchange interaction. It must be responsible for the softening 
of phonons and kinks in the dispersion relation of quasi-particles. However, the superexchange 
interaction must be mainly responsible for the formation of Cooper pairs. 

PACS numbers: 74.20.-z,71.38.-k, 75.30.Et 
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It is an important issue to elucidate the mechanism of 
high-T c superconductivity occurring on CuO? planes [l| . 
Two observations, the softening of phonons \A Li IJ. la, lol 
and kinks in the quasi-particle dispersion pf la, |3, llflj - 
imply the relevance of the electron-phonon interaction. 
One may argue that it must be responsible for high-T c 
superconductivity Its origin and role should be clar- 
ified. 

Doped holes mainly go into O ions. This implies that 
the charge susceptibility of 3c? electrons on Cu ions is 
much smaller than that of 2p electrons on O ions and 
charge fluctuation of 3d electrons can never be devel- 
oped. Then, the conventional electron-phonon interac- 
tion, which directly couples with charge fluctuations, can 
play no crucial role. On the other hand, antiferromag- 
netic spin (AFS) and superconducting (SC) fluctuations 
are certainly developed. One may argue that their devel- 
opments are because of the su pere xchange interaction. It 
is shown in early papers |lll Il2fl that the condensation 
of G?7~wave Cooper pairs bound by the superexchange in- 
teraction can explain observed T c . It is shown in a previ- 
ous paper |lj| that pseudo-gaps appear because of large 
life-time widths of quasi-particles due to SC fluctuations. 
Suggested by these arguments, we propose in this Letter 
an electron-phonon interaction arising from the modula- 
tion of the superexchange interaction by phonons. 

It is shown in another previous paper Q| that 
Gutzwiller's quasi-particle band Il5l lies between the 
lower and upper Hubbard bands [Ifj in metallic phases 
in the vicinity of the Mott-Hubbard transition. The 
superexchange interaction arises from the virtual ex- 
change of pair excitations of electrons in spin channels 
across the lower and upper Hubbard bands; as long as 
the Hubbard splitting is significant, it works between 
Gutzwiller's quasi-particles or their renormalized ones 
|l7j . When nonzero bandwidths of the lower and up- 
per Hubbard bands are ignored, the exchange constant 
between nearest-neighbor Cu ions is given by 



J 



W 4 



{e d + U- e p f 



1 



(1) 



with V the transfer integral between 3d and 2p orbits 
on adjacent Cu and O ions, and e d and e p their energy 



levels. The exchange constant depends on V, e p and e^'s 
of adjacent Cu ions in such a way that 
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with Ae d i a variation of e d of the ith Cu ion, Ae p ryi a 
variation of e p of the [ij]th O ion that lies between the 
ith and jth Cu ions, and AVi m that of V between the 
ith Cu ion and the [ij]th O ion. When we take the x- 
and y-axes along Cu-O-Cu bonds, they are given by 
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with A d , A p and Ay being constants, Ui the displacement 
of the ith Cu ion, u,.£ s that of an O ion on the adjacent 
s = + or s = — side along the £-axis of the ith Cu ion, 
e x = (1,0), e y = (0, 1), and = (R, - R^/IR, - Rj|, 
with Rj the position of the ith Cu ion. Displacements of 
the ith Cu and the [ij]th O ions are given by 
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with R[yi = (Ri+Rj)/2, M d the mass of Cu ions, M p the 
mass of O ions, 6 Aq and 6 A q annihilation and creation 
operators of phonons with polarization A and wave vector 
q, uj\q energies of phonons, e A q unit polarization vectors, 
and -/V the number of unit cells. The q dependence of 
v d ,\q and fp.Aq can play a crucial role; v d ,\q = for 
the breathing modes that bring no changes in adjacent 
Cu-Cu distances while «p,Aq — 1 for such modes. 

The electronic part can be well described by the t- 
J model on a square lattice: TC = — ^2ij a tijd\ a dj a — 
(J/2)J2(ij)( s i ■ s j) + UooJ2i n i] n ih witn tne summa- 
tion over (ij) restricted to nearest neighbors, Sj = 
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Pauli matrixes, and rii a = d\ di a . An infinitely large on- 
site repulsion, i/oo/lfyy) I — > +oo, is introduced to exclude 
any doubly occupied sites. According to Eq. J3J|, there 
are two types of electron-phonon interactions. When 
only longitudinal phonons are considered or when e^q = 
{Qx>Qy>Qz)/q is assumed, they are given by 
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%(q) = 2 [(q x /q) sin (q x a/2) + (q y /q) sin (q y a/2)}, %(q) 
= 2 [(fe/g) sin(g x a/2) - sin (g„a/2)], rfe(k) = 

cos(fc a; a) + cos(fc y a), and %(k) = cos(fc x a) — cos(fc y a), 
with a the lattice constant. 

We follow the previous paper 0] to treat the infinitely 
large Uoo, where a theory of Kondo lattice is developed. 
A renormalized single-site approximation (SSA), which 
includes not only all the single-site terms but also the 
Fock term AE(k) due to the superexchange interaction, 
is reduced to solving the Anderson model. The self- 
energy of the Anderson model is expanded as T, a (ie n ) — 



E(0) + (l-^)ie n +E^(l- 



<)A^ 



small spin-dependent chemical potential shift. Note that 
7 = <$> aa . The Wilson ratio is defined by W s = ffis/^f, 
with (p s = ^> a(J — ^) a _ a . For almost half filling, charge fluc- 
tuations are suppressed so that <fi c — 4> aa + 4> a ~ a <C 1. 
For such filling, ^ > 1 so that <^> s ~ 2^> 7 or W s a 2. 
The dispersion relation of quasi-particles is given by 
£(k) = (l/0 7 ) [- E j l,,' ,k n " + S(0) + AE(k) - /i] , 
with fi the chemical potential. 

The irreducible polarization function ir s (iu>i , q) in spin 
channels is divided into single-site Tr s (iu)i) and multi- 
site Air s (iu>i , q) . The spin susceptibility is given by 
X s (jw;,q) = 2n s (iuji,q)/ {l - [|j(q) + t/oc] 7T s (ia^,q)}, 
with J(q) = 2J?7 s (q). In Kondo lattices, local spin fluc- 
tuations at different sites interact with each other by an 
exchange interaction. Following this physical picture, we 
define an exchange interaction I s (iun , q) by 

Xsfah q) = Xsfai)/ [l - \ls(iui,q)Xs(iui)] , (7) 

with x s (iuji) = 2tt s (ioJi) / [1 — Uao%s(iwi)] the susceptibil- 
ity for the Anderson model. Then, I a (ioji,a} = J(q) + 
2U% c Air s (iuJi, q). When the Ward relation is made 



use of, the irreducible single-site three-point vertex func- 
tion in spin channels, \ s (ie n ,ie n -\-iiOi;iuii), is given by 



lim UocXsfenjiSn + uji;kot) = 2(f> s /x s (iu>i), (8) 

for |e„| — > +0 and |u>;| — > +0. We approximately use 
Eq. © for |e„| < k B T K and < k B T K , with T K the 
Kondo temperature defined by k B Tn = [l/x s (0)] T ^ o . 
The main term of 217^ Air s (iuii, q) is an exchange interac- 
tion arising from the virtual exchange of pair excitations 
of quasi-particles, which is less effective than J(q). The 
so called spin-fluctuation mediated interaction, whose 
single-site term should be subtracted because it is con- 
sidered in SSA, is given by ±[2</> s /x s (iw;)] [x s («w;,q)- 
Xs(iu>i)]. It is simply given by ^i/*(iw/,q), with 
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-7*(iwj,q) = 



^/ s (^,q) 
\l s (iuJi,q)x s (iu>i)' 



(9) 



Because of these equations, we call 7 s (io;j,q) a bare ex- 
change interaction, I*(iu)i,ci) an enhanced one, and <f) s 
an effective three-point vertex function in spin channels. 

The enhanced one is expanded as I* (iwi , q) = Jq + 
2/*77 s (q) + • • • . The nearest-neighbor Jj" is mainly re- 
sponsible for the development of not only SC but also 
charge bond-order (CBO) fluctuations |l9j. Because con- 
tributions from \ui\ < k B Tx are the most effective, we 
ignore its energy dependence. An effective SC suscep- 
tibility, which is multiplied by <p 2 s , is calculated in the 
ladder approximation with respect to Jj" : 
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with Api a a for r = d wave, with 
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Only d-wave SC fluctuation are considered in this Letter. 
An effective CBO susceptibility is similarly given by 
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for r = s, p and ci waves, with 



(CBO)/- \ k B T 2/1 \ 1 



k k„-^(k-|q) 
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The form factors of p waves are defined by f] x (k.) 
y/2sin(k x a) and %(k) = v2sin(A:ya). 



FIG. 1: Four processes contributing to the renormalization of phonons. A solid line stands for an electron, a broken line for a 
phonon, a wavy line for the superexchange interaction J, and a solid circle for the effective vertex function <j> 3 . 



A renormalized Green function for phonons is given by D\(iuji,q) — 2uj\ ci / [(icui) 2 — Lu^ x + 2ui\qAuJx(iuji, q)] , with 
Acu\(iuJi,q) — — (h 2 /2M p w Aq ) S(iuJi,q). Because phonons are renormalized by AFS, SC and CBO fluctuations as 
well as pair excitations of quasi-particles in charge channels or charge density fluctuations, we consider four processes 
shown in Fig.^ When only the part of T = s in Eq. is considered, it follows that S — S s + S sc + S CBO + S c , with 

S s («w ; ,q) = ^Y^(q L )^-J2i 1 2 s (q')x s (iuJi+iuJi l ,ci' + iq) Xsfai ~ iu) V , q' + |q) , (14a) 
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S sc (iw ; ,q) = ^Y^(q)^-J2x S d C (m + ^vM + |q) xTfal - *wj',-q' + |q) , (14b) 
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S CBO (^,q) = J3 Y A 2 (q) Y, ^fr E xP°(^ + iwi'.q' + |q) X? B °(^ - -q' + |q) , (14c) 

■' ' i'q' 
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with 1a (q) = ?/s(q) Cp^Aq^sdq) + CdVd,\q\/ M p /Md ■ Here, Z(ie n , iu>i; k, q) is a vertex function in the charge- 
density channel. It is also enhanced by AFS, SC and CBO fluctuations; Z = Z s + Z sc + Z CBO + • • • , with 

Z a «e n ,*a;,;k,q = Y A (q) — r 2^ r? s q' > . , . ; ~V ^77 / , 1 x ' ( 15a ) 

N ^ i£n + iu[ + ±ui - £(k + q' + ^q) 

Z sc (ie n ,icj J ;k,q) = ^(q)^^^ 7 ^ ( k ~ 5 q ' + 3 q ) % ( k ~ 5 q ' + 3 q ) 

i'q' 



iQ c (iuy - t^wj.q' - |q) #f + i|«i,q' + £q) 



-ie„ + iwj - ifyWi - £(-k + q' - |q) 
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(15b) 



Z CBO (ie n ,iui;k,q) = -r^fq)-^^^* (k + |q' + iq) ?7r (k+ |q' + |q) 



with 



i'q' r 

^ i^ BO + il^.q' + lq) i^ B ° (-t^y + t|t^ q' + |q) 
ie„ + icJ; + iiw; - £(k + q' + iq) 

j, { . x 1 ^sc r x -|/r7rf^,q) -|jr7rg BO (z^,q) 

# s («w;,q) = - — 5— — # d (twi,q)= \ -, K T {m,<V = - \ CBO ,. r- (16) 

No softening occurs for q = because ?7 s (q — > 0) oc |q|. When q goes from T point to the zone boundary, 
the softening must increase first. However, it is unlikely that the softening is the largest at the zone boundary. 
For example, consider the breathing mode at X point, qx = (±7r/a,0). Because Vd,\q x = 0, the electron-phonon 
interaction described by Eq. (|5b|) vanishes. For q ^ q^, Vd,\q is nonzero. This implies that the softening may not be 
the largest at X point alongT-X line. In actual, several experimental data show that the softening is the largest for 
q a little different from qx 0, E| . 

Because a low-energy scale is ksTx, we put 

S(iu>i,q) ~ s/ksTK, (17) with s — 0(1) a dimensionless constant. It follows that 

Auj x (uJ\q x ,(lx) ^ -15sc 2 meV, (18) 



4 



for 



C p ~ c p eV/A, w Aqx = 0.1 eV, k B T K = 0.1 eV, (19) 
with c p = O(l) a dimensionless constant. When we take 



V ~ 1.6 eV, e 3rf - e 2p ^ -1 eV, C/ ~ 5 eV, 



(20) 



following the previous paper (l7|. it follows that c p = 
(0.5-1) for A d = (1-2) eV/A. If sc 2 p ~ 1, the observed 
softening as large as 0- C3- [3- 0] 



Aw A (wA qx ,qx) - -10 meV, 



(21) 



can be explained. It should be examined whether sc p is 
actually as large as 1. 

A process corresponding to Fig. 1(d) renormalizcs 
quasi-particles. The self-energy correction is given by 

J-AE(ie„,k) = --f^ VV( ie n, iwf, k, q)D\(iui, q) 

<^ ^ AZq 

x- . 1 ... r, (22) 



ie n + iui - £(k + q) ' 



with 



g\(ie n ,i(jJi;\t,q) = C p 
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W^Z(ie n , iwi;k,q). 



(23) 

It is likely that the contribution of Fig. [Jd) dominate 
those of the other three, Figs. Qla)-(c). I n such a case, 



than fl/fl 



100 meV. The main Cooper-pair interac- 



tion must be the superexchange interaction. 

There are two other types of electron-phonon interac- 
tions. Note that (f> c and l/</> 7 are small parameters in 
the vicinity of the Mott-Hubbard transition. The con- 
ventional one arising from the modulation of 3d-electron 
levels, which can directly couples with charge fluctua- 
tions, gives renormalization effects higher order in <p c and 
1/^7 . The interaction arising from the modulation of Uj 
gives renormalization effects higher order in l/0 7 . Then, 
we ignore both of them. On the other hand, what are 
considered in this Letter are of the order of (4> c )° (1 / 4>-y)° • 

In conclusion, we propose the electron-phonon interac- 
tion arising from the modulation of the superexchange in- 
teraction by phonons, which is only relevant for strongly 
correlated electron liquids in the vicinity of the Mott- 
Hubbard transition. Its novel property is that it can be 
enhanced by spin, superconducting, and charge bond- 
order fluctuations as well as charge density fluctuations. 
A phenomenological argument where parameters are de- 
termined from the observed softening of phonons implies 
that the enhanced electron-phonon interaction is also re- 
sponsible for kinks in the dispersion relation of quasi- 
particles in cuprate-oxide high-T c superconductors. How- 
ever, it can never be the main Cooper-pair interaction. 
The main one must be the superexchange interaction. 

This work is supported by a Grant-in- Aid for Scientific 
Research (C) No. 13640342 from the Japan Society for 
the promotion of Science. 



g\(ie n , iw/;k, q) ~ J k B T K \ALJx(^x<i x > <lx)l- ( 24 ) 

When the experimental value (|21|l is used, we obtain 

g\{ie n , icui;k, q) ~ 30 meV. (25) 

This is large enough for optical phonons to cause kinks 
in the dispersion relation of quasi-particles. Two types of 
kinks are observed fl(3 | . The renormalization by phonons 
can explain one type of kinks observed in both normal 
and SC phases. However, it cannot explain the other 
type of kinks observed only in SC phases. 

The phonon-mediated pair interaction is given by 
ff 2(0,0;k,q)^A(0,q) or -2 5 j[(0,0;k > q)/w Aq . The soft- 
ening of phonons is the largest for q ~ qx along T-X 
line. This implies that the pair interaction by phonons 
is attractive between nearest neighbors. The nearest- 
neighbor part of — 2<7a(0, 0; k, q)/wAq should be included 
in I J*. According to the argument in this Letter, it fol- 
lows that 



-2^(0,0;k,q x )/ W Aq x =s -20 meV. 



(26) 



The phonon-mediated interaction cannot be ignored in 
cuprate-oxide superconductors. However, it is smaller 
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